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I.  INTRODUCTION 


The  transonic  flow  over  bodies  of  revolution  is  of  considerable  interest 
in  aerodynamics.  The  critical  aerodynamic  behavior  occurs  in  this  speed 
regime  and  a  rapid  increase  in  the  aerodynamic  coefficients  is  observed.  The 
flowfield  is  complex  due  to  the  presence  of  shocks  and  can  be  characterized  by 
strong  vi scous-i nvi sci d,  shock-boundary  layer  interactions.  Further  complica¬ 
tions  often  occur  in  practical  flows  of  interest  due  to  the  presence  of  sepa¬ 
rated  flow  regions  which  is  especially  true  for  flow  at  moderate  to  high 
angles  of  attack.  As  a  result,  it  is  desirable  to  use  the  Navier-Stokes  comp¬ 
utational  technique  to  compute  such  flows. 

In  recent  years  Navier-Stokes1" 3  computational  methods  have  been  used  to 
compute  flow  over  bodies  of  revolution  at  transonic  speeds.  In  Reference  1 
the  three  dimensional  flow  over  a  boattailed  afterbody  was  computed  where  the 
primary  emphasis  was  on  the  boattail  flow  field.  Calculations  of  three  dimen¬ 
sional  flow  over  a  projectile  have  been  made  in  Reference  2;  however,  the  wake 
or  base  region  flow  was  not  computed.  The  wake  region  has  been  included  in  a 
three  dimensional  flowfield  computation  in  Reference  3.  These  calculations 
were  made  with  the  compressible,  thin-layer  Navier-Stokes  equations  which  were 
solved  using  the  implicit  Beam  and  Warming  central  finite  difference 
scheme.4”6  Generally,  these  calculations  did  not  have  sufficient  grid  resolu¬ 
tion  due  to  lack  of  adequate  computer  resources.  With  the  advent  of  computers 
such  as  the  CRAY  X-MP/48  and  CRAY  2  and  the  availability  of  the  large  memory 
(core  memory  of  CRAY  2  and  fast  external  memory  device,  SSD  of  the  CRAY  X- 
MP/48),  it  is  now  possible  to  use  fine  computational  grids  needed  for  three 
dimensional  numerical  simulation  of  transonic  flows.  Reference  7  is  an 
excellent  example  where  the  CRAY  2  computer  was  used  to  perform  transonic  flow 
computations  over  a  hemisphere  cylinder  at  high  angles  of  attack  using  its 
large  incore  memory.  In  the  present  work,  we  make  extensive  use  of  the  SSD 
device  of  the  CRAY  X-MP/48  to  provide  sufficient  grid  resolution.  This  neces¬ 
sitates  splitting  the  data  base  of  a  large  single  grid  into  smaller  pieces 
that  fit  within  the  incore  memory  of  the  X-MP/48.  While  this  break  up  can  be 
achieved  in  various  ways,  the  simplest  approach  is  to  break  the  single  grid 
into  a  number  of  smaller  grids. 

The  numerical  scheme  used  is  an  implicit  scheme  based  on  flux-splitting8 
and  upwind  spatial  differencing  in  the  streamwise  direction.  Upwind  schemes 
can  have  several  advantages  over  central  difference  schemes,  including  natural 
numerical  dissipation  and  better  stability  properties. 

The  accuracy  of  the  composite  grid  scheme  has  been  tested  by  numerically 
simulating  the  transonic  flow  over  an  ellipsoid  at  10°  angle  of  attack  and 
comparing  this  result  with  the  result  obtained  with  a  single  grid.  In  addi¬ 
tion,  the  flow  over  a  secant-ogi ve-cyl inder-boattai 1  projectile  has  been 
computed  using  a  fine  grid  at  M  =  D.gfj,  l.i  and  4°  angle  of  attack.  These 
computed  results  have  been  compared  with  experimental  data  to  determine  the 
accuracy  of  the  numerical  predictions.  In  Section  II,  the  governing  equa¬ 
tions,  numerical  procedure  as  well  as  the  composite  grid  scheme  have  been 
described.  Results  are  presented  in  Section  III. 


II.  NUMERICAL  METHOD 


1.  GOVERNING  EQUATIONS 

The  three  dimensional  Navier-Stokes  conservation  equations  of  mass, 
momentum,  and  energy  can  be  represented  in  flux  vector  form  as: 

+  VF'  +  Fv>  +  VG  +  Gv)  +  V"  +  Hv}  =  0  (1) 


where  the  independent  variable  t  is  the  time  and  the  spatial  variables  C,n,C 
are  chosen  to  map  a  curvilinear  body  conforming  discretization  into  a  uniform 

computational  space.  Here  Q  contains  all  the  dependent  variables  and  F,  G  and 

A  A  A  A 

H  are  the  inviscid  fluxes.  The  flux  terms  Fv,Gv  and  Hy  contain  viscous  deriv¬ 
atives  and  throughout  a  nondimensi onal  form  of  the  equations  is  used.  The 
conservative  form  of  the  equations  is  maintained  chiefly  to  capture  the 
Rankine  Hugoniot  shock  jump  relations  as  accurately  as  possible. 

For  body  conforming  coordinates  and  high  Reynolds  number  flow,  if  c  is  the 
coordinate  away  from  the  surface,  the  thin  layer  approximation  can  be  made  in 
the  c  direction  and  the  governing  equations  can  be  written  as: 


3  Q  +  3  F  +  3  G  +  3  H  =  Re_13_S  .  (2) 

^  s  n  £  s 

Here  the  viscous  terms  in  c  have  been  collected  into  the  vector  S  and  the 
nondimensional  reciprocal  Reynolds  number  is  extracted  to  indicate  a  viscous 
flux  term. 

In  differencing  these  equations  it  is  often  advantageous  to  difference 
about  a  known  base  solution  denoted  by  subscript  o  as: 


-1 


-  V  +  VF  -  +  VG  -  Go)  +  VH  +  V  -  Re  VS  "  V 


=  -3  Q  -  3.F-  -  3  G  -  3  H  +  Re_i3  S 
t'o  CO  no  co  CO 


Ci) 


where  6  indicates  a  general  difference  operator,  and  3  is  the  differential 
operator.  If  the  base  state  can  be  properly  chosen,  the  differenced  quanti¬ 
ties  can  have  smaller  and  smoother  variation  and  therefore  less  differencing 
error.  The  freestream  is  used  as  a  base  solution  in  the  present  formulation. 


2.  IMPLICIT  FINITE  DIFFERENCE  ALGORITHM 

The  implicit  approximately  factored  scheme  for  the  thin  layer  Navier- 
Stokes  equations  that  uses  central  differencing  in  the  n  and  c  directions  and 
upwinding  in  C  is  written  in  the  form: 


* 
♦ 


[I  +  h<sy(A+)n  +  h<$  Cn  -  hRe'1_S  J_1MnJ  -  D .  |  ] 

S  4  4  Is 


d  +  h6j( A‘)n  +  h6nBn  -  n.|n]AQn  = 


-At{5^|(F+)n  -  F+]  +  6j|(F')n  -  F;]  +  6n(Gn  -  GJ 


+  6c(Hn  -  Hj  -  Re-%;(Sn  -  SJ>  -  De(Qn  -  6.) 


where  h  =  At  or  (At)/2  and  the  freestream  base  solution  is  used.  Here  6  is 


typically  a  three  point  second  order  accurate  central  difference  operator,  6 

is  the  midpoint  operator  used  with  the  viscous  terms,  and  the  operators  6^ 

f  ^  ~ 

and  are  backward  and  forward  three-point  difference  operators.  The  flux  F 

has  been  split  into  F  and  F  ,  according  to  its  eigenvalues  and  the  matrices 


A,B,C  and  M  result  from  local  linearization  of  the  fluxes  about  the  previous 
time  level.  Here  J  denotes  the  Jacobian  of  the  coordinate  transformation. 
Dissipation  operators,  ne  and  D-j  are  used  in  the  central  space  differencing 
directions. 7 


The  factored  left  hand  side  operators  can  be  readily  inverted  by  sweeping 
and  inversion  of  tridiagonal  matrices  with  5  *  5  blocks.  This  two  factor 
implicit  scheme  is  readily  vectorized  or  multi-tasked  in  planes  of  S  = 
constant. 


3.  COMPOSITE  GRID  SCHEME 


In  the  present  work,  a  composite  grid  scheme  has  been  developed  in  which  a 
large  single  grid  is  split  into  a  number  of  smaller  grids  so  that  each  of 
these  grids  can  use,  in  turn,  the  available  core  memory  while  the  rest  are 
stored  on  the  SSD  device.  The  use  of  such  composite  scheme  however,  requires 
special  care  in  storing  and  fetching  the  interface  boundary  data. 


If  a  code  is  well-written  in  a  modular  form  for  a  single  grid,  it  is  a 
relatively  simple  task  to  reconfigure  it  either  as  a  composite  grid  code,  or  a 
blocked  grid  data  structure  code  using  an  external  large  memory  device  such  as 
the  SSD  of  the  CRAY  X-MP/48.  The  simplest  case  to  begin  with  is  the  blocked 
grid  data  structure  code  in  which  a  single  large  grid  is  partitioned  into  a 
series  of  smaller  grid  blocks  as  sketched  in  Figure  1.  Such  a  partitioning  of 
a  single  original  grid  may  come  about  because  all  of  the  single  grid  will  not 
fit  into  a  particular  computer's  high  speed  memory  or  because  a  block  itera¬ 
tive  technique  is  in  use.  To  make  a  code  work  for  a  partitioned  grid  that 
works  for  a  single  grid  it  is  necessary  to  build  an  external  data  management 
scheme  for  a  block  update  process  (e.g.  block  iterative  or  block  time  depend¬ 
ent  process).  In  a  block  update  scheme,  if  data  from  one  subgrid  is  used  to 
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update  the  interface  boundary  of  another  subgrid,  the  subblock  solution  will 
approach  the  solution  of  the  single  grid.  (We  say  approach  because  the 
differencing  approximations  near  a  boundary  may  drop  in  order  of  accuracy  or 
use  metrics  evaluated  in  a  different  manner.  Likewise,  iterative  convergence 
may  be  less  tight  in  a  given  block.) 

A  code  developed  for  a  single  grid  can  be  made  to  work  for  a  block  grid 
structure  by:  1)  mapping  and  storing  the  information  for  each  grid  onto  a 
large  memory;  and  2)  supplying  interface  boundary  arrays,  pointers  and 
updating  procedures.  Consider  the  situation  in  Figure  1  in  which  the  single 
grid  from  J  =  1,  Jmax  is  partitioned  into  three  grids.  Although  not  shown  in 

this  figure,  the  base  region  of  the  projectile  is  easily  included  by  adding 

one  or  more  zones  in  that  region.  We  begin  the  procedure  with  grid  Gl.  For 

each  grid,  an  estimated  solution  (or  initial  quess)  is  put  into  the  large 

memory  along  with  all  pertinent  information  such  as  the  grid,  metrics,  M  ,  R 

00  0 

etc.  This  data  is  brought  into  the  working  memory  and  appropriate  boundary 
update  information  is  also  called  in  and  used  to  update  any  interface  bound¬ 
aries.  At  this  point,  the  problem  of  dealing  with  a  single  grid  code  must  be 
dealt  with.  Solution  variables  are  updated  as  before  (more  than  one  iterative 
update  per  subgrid  is  possible  or  interface  boundaries  can  be  iteratively 
corrected).  The  updated  solution  is  then  stored  back  into  large  (external) 
memory.  Also  stored  into  the  appropriate  locations  is  any  information  needed 
to  update  other  interface  boundaries.  This  procedure  is  repeated  in  the  next 
subgrid  and  so  on. 

The  chief  complexity  of  any  blocked  or  composite  grid  scheme  is  the  pro¬ 
blem  of  storing  and  fetching  the  interface  boundary  data  from  its  appropriate 
place,  i.e.,  settinq  up  the  proper  bookkeeping.  Additional  complications  come 
about  when  we  try  to  make  blocked  grids  with  interface  boundaries  give  a  solu¬ 
tion  which  is  identical  to  the  solution  obtained  by  a  large  single  grid  with 
no  interface  boundaries.  This  involves  maintaining  the  same  or  equivalently 
good  metrics,  differencing  accuracy,  and  conservation  statements  etc.,  at  or 
near  interface  boundaries  which  would  exist  on  a  continuous  smooth  grid.  One 
choice  that  allowed  a  simpler  management  of  the  interface  boundary  data  while 
still  maintaining  second  order  accuracy  was  to  use  central  finite  differences 
at  these  boundaries.  As  will  be  shown  in  the  results,  this  procedure  works 
rather  well  compared  to  dropping  the  accuracy  of  the  upwind  scheme  from  second 
order  to  first  order. 

For  the  simple  partitioning  shown  in  Figure  1,  all  subgrid  points  are 
members  of  the  original  grid  and  no  interpolations  are  required.  This  pro¬ 
cedure  thus,  has  the  advantage  over  patched  or  overset  grid  schemes  which  do 
require  interpolations.  The  partitioned  grid  has  four  interface  boundaries, 
Jl  =  Jimax»  ^2  =  J2  =  J2max  and  ^3  =  Data  f°r  these  planes  are  to  be 

supplied  from  the  other  grids  by  injecting  interior  values  of  the  other  grid 
onto  the  interface  boundaries.  For  example,  the  =  Jimax  boundary  of  Gl  may 

be  taken  as  the  J2  =  2  or  the  =  3  plane  interior  values  of  G2.  More  over¬ 
lap  would  mean  more  implicit  like  behavior.  Because  only  one  grid  and  some 
boundary  and  pointer  arrays  can  be  in  the  working  memory  at  one  time,  a  low 


storage  scheme  to  shuffle  data  from  one  grid  to  another  must  be  used.  In 
working  on  grid  Gl  all  of  the  usual  data  Q,x,y,z...  of  grid  Gl  plus  the 
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interface  boundary  data  stored  in  an  array,  say,  OB  from  grid  G2,G3,...  will 
have  to  be  in  working  memory.  The  array  QB  need  only  store  Q  on  the  interface 
boundary  as  other  data  such  as  metrics  on  the  interface  boundary  can  be  compu¬ 
ted  from  the  usual  data  x,y,z...  After  updating  grid  Gl,  interface  boundary 
data  QB  supplied  by  Gl  for  the  other  grids  should  be  loaded  while  all  data  of 
Gl  is  still  in  memory.  Thus  the  data  array  QB  should  be  shuffled  back  into  a 
permanent  large  external  memory,  say  QBC,  in  such  a  way  that  it  can  be  readily 
brought  back  for  grid  G2,  grid  G3,  etc.  There  are  several  ways  to  store  QB  in 
memory,  the  approach  selected  here  is  flexible  and  is  used  in  the  Chimera 
code. 9 

III.  RESULTS 

The  implicit  time  marching  procedure  was  used  to  obtain  the  desired  steady 
state  result  starting  from  initial  freestream  conditions  everywhere.  Boundary 
conditions  were  updated  explicitly  at  each  time  step.  The  solution  residual 
dropped  at  least  three  orders  of  magnitude  before  converged  solutions  were 
obtained.  In  addition,  the  surface  pressure  distribution  was  checked  for  time 
invariance.  Solutions  have  been  obtained  for  two  cases:  (i)  flow  over  an 
ellipsoid  at  =  0.90  and  a  =10°,  primarily  to  check  the  accuracy  of  the  com¬ 


posite  grid  solution  and  (ii)  flow  over  a  projectile  at  M0 
4°  which  is  calculated  using  a  fine  grid. 


0.96,  1.1  and  a  = 


First,  the  flow  over  an  ellipsoid  at  =  0.90  and  10°  angle  of  attack  is 

considered.  Figure  2  shows  the  computational  grid  for  the  ellipsoid.  Figure 
2a  shows  the  surface  of  the  ellipsoid  while  Figure  2b  is  a  cross-section  of 
the  full  grid  showing  the  longitudinal  grid  point  distribution.  This  grid  has 
41  points  in  the  streamwise  direction,  IB  points  in  the  circumferential,  and 
31  points  in  the  normal  direction.  Such  a  small  grid  was  chosen  so  that  the 
full  grid  solution  could  be  obtained  using  a  CRAY  X-MP/48  with  a  two  million 
word  central  core  memory.  The  full  grid  was  then  partitioned  into  three 
smaller  grids  and  the  SSH  device  was  used  to  obtain  the  solution  for  this 
composite  grid  scheme.  Figure  3  shows  the  surface  pressure  distribution  in 
the  lee  side.  Here,  the  full  grid  solution  is  shown  as  a  solid  line  and  the 
composite  grid  solution  as  a  dotted  line.  The  differences  in  these  curves  are 
apparent  and  occur  near  X  =  .3  and  .8  at  the  interfaces  between  the  parti¬ 
tioned  grids.  The  spatial  accuracy  at  the  interface  was  first  order  in  this 
case.  Improving  the  spatial  differencing  accuracy  at  these  boundaries  to 
second  order  resulted  in  the  solution  shown  in  Figure  4.  As  shown  in  this 
figure,  the  composite  grid  solution  is  virtually  identical  to  the  full  grid 
solution. 


Second,  the  flow  over  a  projectile  at  a  transonic  speed  =  0.96  and  at  a 


higher  transonic  speed  M  =  1.1  for  a  = 


was  computed.  The  model  used  for 


the  experiment  and  computational  study  presented  here  is  an  idealization  of  a 
realistic  artillery  projectile  geometry.  The  experimental  model  shown  in 
Figure  5  consists  of  a  three-cal i ber  (one-caliber  =  maximum  body  diameter), 
sharp,  secant-ogive  nose,  a  two-caliber  cylindrical  mid-section  and  a  one- 
caliber  7°  conical  afterbody  or  boattail.  A  similar  model  was  used  for  the 
computational  studies  with  the  only  difference  being  a  6°;  rounding  of  the  nose 
tip.  The  nose  tip  rounding  was  done  for  computational  efficiency  and  is 
considered  to  have  little  impact  on  the  final  integrated  forces. 


Experimental  pressure  data10’11  are  available  for  this  shape  and  were 
obtained  in  the  NASA  Langley  eight  foot  Pressure  Tunnel  using  a  sting  mounted 
model.  The  test  conditions  of  1  atm  supply  pressure  and  320  K  supply  tempera¬ 
ture  resulted  in  a  Reynolds  number  of  4.5  x  io6  based  on  model  length. 

The  computational  grid  used  for  this  computation  is  shown  in  Figure  6. 
Figure  6a  shows  the  longitudinal  cross  section  of  the  3D  grid  while  Figure  6b 
shows  the  ci rcumferenti al  distribution  of  grid  points.  As  shown  in  Figure  6a, 
the  clustering  of  grid  points  near  the  body  surface  is  done  to  resolve  the 
viscous  boundary  layer  near  the  body  surface.  Grid  clustering  has  also  been 
used  in  the  longitudinal  direction  near  the  boattail  and  the  base  corners 
where  large  gradients  in  the  flow  variables  are  expected.  The  grid  consists 
of  162  points  in  the  streamwise  direction,  36  points  in  the  circumferential 
direction  and  50  points  in  the  normal  direction.  This  amounts  to  a  little 
over  12  million  words  of  storage  for  the  code  on  the  CRAY  X-MP/48.  Only  up  to 
4  Mw  of  central  core  memory  was  easily  accessible;  therefore,  the  full  grid 
was  partitioned  into  four  smaller  grids  each  of  which  would  use  the  core 
memory  in  turn  while  the  rest  is  stored  on  the  SSD  device.  In  addition,  the 
grid  is  made  to  wrap  around  the  base  corner  to  facilitate  the  generation  of 
the  large  single  grid.  What  affect  it  has  on  the  base  region  flow  is  current¬ 
ly  under  investigation  and  has  not  been  studied  here.  The  primary  emphasis  in 
this  research  was  on  the  use  of  the  composite  grid  technique  itself.  The 
computations  were  performed  on  the  CRAY  X-MP/48  at  the  US  Army  Ballistic 
Research  Laboratory. 

For  the  computation  of  turbulent  flow,  a  turbulence  model  must  be 
supplied.  In  the  present  calculation,  a  two  layer  algebraic  eddy  viscosity 
model  due  to  Baldwin  and  Lomax12  is  used.  Figure  7  and  8  show  the  pressure 
contours  and  Mach  contours  for  the  projectile  in  the  windward  and  leeward 
planes.  As  shown  in  Figure  7,  the  pressure  contours  show  the  expansions  at 
the  ogive-cylinder,  boattail  and  base  corners.  The  expansions  at  the  base 
corners  are  rather  wear,.  This  figure  also  indicates  the  presence  of  shock 
waves  on  the  cylinder  and  also  on  the  boattail.  Figure  8  shows  the  Mach 
number  contours.  It  clearly  shows  the  expansions  and  shock  structure  that 
typically  occur  on  the  projectile  at  transonic  speeds.  Sharp  shocks  are 
clearly  observed  on  the  boattailed  flow  field  which  are  also  asymmetrically 
located  (the  one  on  the  wind  side  being  closer  to  the  base  than  its  counter¬ 
part  on  the  lee  side.)  This  figure  also  indicates  the  presence  of  strong  free 
shear  layers  in  the  base  region  flow  field  and  also  shows  the  expected  asym¬ 
metric  flow  pattern  in  the  wake.  The  asymmetric  wake  flow  can  be  better  seen 
in  Figure  8  which  shows  the  velocity  field  in  the  near  wake  for  both  the  wind 
as  well  as  the  leeside.  The  reci rcul atory  flow  pattern  in  the  base  region  is 
evident.  Three  pairs  of  separated  flow  bubbles  can  be  clearly  seen  in  this 
figure.  Each  pair  consists  of  two  vortices,  one  clockwise  and  one  counter¬ 
clockwise.  Two  of  these  pairs  occur  in  the  leeside  whereas  only  one  pair  can 
be  seen  in  the  windside.  The  experimental  model  was  sting  mounted  in  the  base 
region  and  current  work  is  in  progress  to  model  the  sting  before  accuracy  of 
the  base  region  flow  prediction  can  be  determined.  In  the  present  computa¬ 
tion,  the  sting  has  not  been  modeled. 

Figure  10  shows  the  surface  pressure  distribution  as  a  function  of  the 
longitudinal  position  for  the  wind  side.  The  computed  result  is  compared  with 
experimental  data10  shown  in  circles  and  is  in  very  good  agreement.  The 
expansions  and  recompressions  near  the  ogive-cylinder  and  cyl i nder-boattai 1 
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junctions  are  captured  adequately  by  the  computations.  Comparison  of  the  sur¬ 
face  pressure  on  the  lee  side  is  shown  in  Figure  11.  A  small  discrepancy  can 
be  seen  on  the  cylinder  and  boattail  sections.  The  computed  result,  however, 
is  generally  in  good  agreement  with  the  experimental  data.  Figures  12  and  13 
show  the  ci rcumferenti al  surface  pressure  distributions  for  longitudinal 
positions  X/D  =  6.19  and  6.56,  respectively.  Both  of  these  longitudinal 
positions  are  on  the  boattail.  In  these  figures,  the  present  computed  result 
is  compared  with  experimental  data  shown  in  circles  and  also  with  a  previous 
computational  result  obtained  on  a  rather  coarse  mesh. 3  Apart  from  the  small 
discrepancy  near  the  lee  side  ($>  =  180°),  the  present  computed  result  is  in 
qood  agreement  with  experimental  data  and  is  more  accurate  compared  to  the 
previous  compiled  result  (see  Figure  12).  The  agreement  of  the  present  result 
w’th  experiment  is  very  good  near  the  wind  side  ($  =  0).  Figure  13  shows  the 
circur#e’-pr*’s’  surface  pressure  comparison  at  a  longitudinal  position  X/D  = 
6.66.  tnis  figure,  the  present  numerical  result  is  in  excellent 

agreemer*  «  * r  »ne  experimental  data.  The  improved  accuracy  of  the  present 
numeri;a'  .  e":*-or'  compared  with  the  previous  coarse  grid  solution  is 
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Nu^erica1  solution  has  also  been  obtained  at  a  high  transonic  speed  Hw  = 

1.1  and  a  =  4- .  Comparisons  of  the  longitudinal  surface  pressure  distribu¬ 
tions  for  Dot h  wind  side  and  lee  side  have  been  made  with  experimental  data10 
and  are  shown  in  Figures  14  and  15.  The  agreement  of  the  computed  surface 
pressures  with  experiment  is  very  good.  At  this  high  transonic  Mach  number, 
the  shocks  on  the  cylinder  as  well  as  on  the  boattail  are  very  weak  as 
evidenced  by  the  absence  of  sharp  rise  in  pressure  in  those  areas.  Strong 
shock  waves  occur  at  lower  transonic  Mach  numbers,  for  example  at  =  0.96 

and  a  =  4°  (see  Figure  10).  The  expansions  and  recompressi ons  near  the  ogive- 
cylinder  and  cyl i nder-boattai 1  junctions  can  be  clearly  observed  in  Figures  14 
and  16. 


IV.  CONCLUDING  REMARKS 

In  conjunction  with  a  new  Navier-Stokes  code,  a  simple  composite  grid 
scheme  has  been  developed  which  allows  fine  computational  grids  needed  for 
accurate  transonic  flow  computations  to  be  obtained  on  a  CRAY  X-MP/48 
computer.  The  numerical  method  uses  an  implicit,  approximately  factored, 
partially  upwind  (flux-split)  algorithm. 

The  accuracy  and  stability  of  the  composite  grid  scheme  have  been  tested 
by  numerically  simulating  the  flow  over  an  ellipsoid  at  angle  of  attack  (M  = 
0.90,  a  =  10°)  and  comparing  this  solution  with  a  single  grid  solution.  The 
flow  fields  over  a  projectile  at  M^  =  0.96,  1.1  and  a  =  4°  have  been  computed 

using  a  fine  grid  and  compared  to  experiment.  The  computed  surface  pressures 
are  in  very  good  agreement  with  experimental  data.  Currently  work  is  in  pro¬ 
gress  to  analyze  and  determine  the  accuracy  of  the  wake  region  flow  and  also 
to  determine  the  critical  aerodynamic  behavior  at  transonic  speeds. 
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Figure  2b.  Longitudinal  cross  section  of  the  computational  grid 
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